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ABSTRACT 

The observational behavior of spherically symmetric inhomogeneous cosmological 
models is studied, which consist of inner and outer homogeneous regions connected by 
a shell or an intermediate self-similar region. It is assumed that the present matter 
density parameter in the inner region is smaller than that in the outer region, and the 
present Hubble parameter in the inner region is larger than that in the outer region. 
Then galaxies in the inner void-like region can be seen to have a bulk motion relative 
to matter in the outer region, when we observe them at a point O deviated from the 
center C of the inner region. Their velocity v p in the CD direction is equal to the 
difference of two Hubble parameters multiplied by the distance between C and O. It is 
found also that the velocity Vd corresponding to CMB dipole anisotropy observed at O 
is by a factor ~ 10 small compared with v p . This behavior of Vd and v p may explain 
the puzzling situation of the cosmic flow of cluster galaxies, when the radius of the 
inner region and the distance CD are about 200 /i -1 Mpc and 40 /i -1 Mpc, respectively 
(Hq = 100/i -1 km sec -1 Mpc -1 ), and when the gaps of density and Hubble parameters 
are « 0.5 and 18%, respectively. 

Subject headings: cosmic microwave background - cosmology: large scale structure of 
the universe - observations 



1. Introduction 

The dipole moment in the cosmic background radiation (CMB) is thought to come mainly from 
the Doppler shift due to the motion of the Local Group (LG), relative to the cosmic homogeneous 
expansion. As the main gravitational source which brings the velocity vector of LG, the existence 
of the Great Attractor (GA) was found by Lynden-Bell et al. (1988) and Dressier (1987). It has 
the position at the redshift of 4300 km sec -1 . On the other hand, the motion of LG in the inertial 
frame consisting of many clusters on larger-scales was studied observationally by several groups: A 
bulk flow of ~ 700 km sec -1 was found by Lauer and Postman (1994) and Colless (1995) as the 
motion of the Abell cluster inertial frame relative to LG in the region with redshift < 15000 km 
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sec -1 , but in the other approach the different result was derived by Giovanelli et al. (1998), Dale 
et al. (1999) and Riess et al. (1997) in the regions with similar redshifts. The Lauer and Postman's 
work is based on the assumption that the brightest cluster galaxies as standard candles and the 
Hoessel relation can be used, but at present these assumptions have been regarded as questionable 
or unreliable. 

Independently of these works, the motion of cluster frames relative to CMB was measured by 
Hudson et al. (1999) and Willick (1999) due to the global Hubble formula using the Tuny-Fisher 
distances of clusters and their redshifts with respect to CMB, and the flow velocity vector was 
derived in the region with about \50h~ 1 Mpc (Hq = 100/t -1 km sec -1 Mpc -1 ). The remarkable 
and puzzling properties of these flows are that the flow velocity reaches a large value ~ 700 km/sec 
on large scale, while the dipole velocity (not due to GA) determined in the form of CMB dipole 
anisotropy seems to be much small, compared with the above flow velocity. 

If the observed large-scale matter motion is caused by the attraction from an over-density 
region containing superclusters, the corresponding velocity must be so large as the large-scale flow 
velocity and it must be reflected in the form of CMB dipole anisotropy. If this motion is caused in 
the spherical void-like region, however, the situation is different, because CMB dipole anisotropy 
near the center can be relatively small in spite of the large-scale flow, as was shown in our previous 
paper (Tomita 1996). In this previous paper an inhomogeneous model on super- horizon scale was 
considered to explain the number evolution of QSOs (Tomita 1995), but the relative smallness of 
the dipole anisotropy can be found independently of the scale of inhomogeneities. The local void 
region was studied independently by Zehavi et al. (1998) as the local Hubble bubble, which has 
the scale ~ 70/i -1 Mpc and is bordered by the dense walls. They analyzed the statistical relation 
between the distances and the local Hubble constants derived from the data of SNIa, and found 
the existence of a void region with a local Hubble constant larger than the global Hubble constant. 
The relation to the SNIa data on larger scales will be discussed from our standpoint in a subsequent 
paper. 

In the present paper we consider more realistic inhomogeneous models on sub-horizon scale, 
corresponding to matter flows ~ 150/j" 1 Mpc, which may be associated with large-scale structures 
or excess powers observed by Broadhurst et al. (1990), Landy et al. (1996), and Einasto et al. 
(1997). In §2, we treat a spherically symmetric inhomogeneous model which consists of inner and 
outer homogeneous regions connected by a shell being a singular layer, and study the behavior of 
large-scale motions caused in the inner region, where the present inner density parameter is smaller 
than the present outer density parameter and the present Hubble parameter in the inner region is 
larger than that in the outer region (a bulk motion in the void-like region was discussed also by 
Nakao et al. (1995)). In this section we treat the single-shell case. The double-shell case and a 
model with an intermediate self-similar region are treated in Appendixes A and B, respectively. In 
§3, we consider light rays which are emitted at the last scattering surface and reach an observer 
situated at a point O deviated from the center C, and CMB dipole and quadrupole anisotropics 
are analyzed. The peculiar velocity of the above large-scale motions and the velocity corresponding 
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to the CMB dipole anisotropy are compared. In §4, a naive explanation about why CMB dipole 
anisotropy is small around the center C is shown. In §5 the consistency of the present models with 
several recent observations of bulk flows is discussed, and in §6 concluding remarks are presented. 



2. Inhomogeneous models and the bulk motion 

In this section we consider spherically symmetric inhomogeneous cosmological models which 
have two homogeneous regions connected with a spherical shell, as shown in Fig. 1. 

The line-elements in the inner region V 1 and the outer region V 11 are described as 

ds 2 = g^dxindxiy = -c 2 (dti) 2 + [a j (t j )} 2 {d( X j ) 2 + [W)] W}, (1) 

where j (= I or II) represents the regions, f 3 ix 3 ) = sm X' 7 )X J an d sinhx- 7 for W = 1,0,-1, 
respectively, and dft 2 = d9 2 + sin 2 9ip 2 . The shell is a time-like hypersurface X given as \ l = Xi 
and x 11 = Xi 1 - 



2.1. Cosmological models 

The Einstein equations are divided into the equations in the two regions and the jump con- 
ditions at the shell. The general formulation of the jump condition at the singular surface was 
derived by Israel (1966) and the concrete expressions of conditions were derived by Maeda (1986) 
and Sakai et al. (1993). Here the expressions by Sakai et al. are shown using the circumferential 
radius of the shell R, the velocity of the shell v 3 , the Lorentz factor j 3 and the Hubble expansion 
parameter H 3 in V (j = I and II) defined by 



R S «'/' = a"/", V S aj %, y> S 1/^1 - Q 2 and W S ^* . (2) 
The Einstein equations for the pressureless matter in the two regions are 

{ Wf + Vc 2 /{ai) 2 = ^ + \Kc 2 , (3) 

where fP is the mass density of matter (oc l/(a J ) 3 ). 

The equations for the surface density a and the velocity v 11 of the shell are expressed as 

~f ll d(4irR 2 a)/dt 11 = UirR 2 -y 2 vp] "'^ (4) 

d(j u v u )/dt u = —j n v u H u + 2-kGo - [ 1 2 v 2 p] u /a, (5) 
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where [tJ?] 11,1 = $ n — The conditions of continuity of the metric (cir 2 = —ds 2 ) and the common 
velocity dR/dr reduce to 

dti/dt 11 = 7V7 11 (6) 

and 

Vt/V + H l R) = 7 u (f' u v u + H ll R), (7) 

where p = df J (x 3 )/dx 3 ■ 

Another important component of jump conditions playing a role of a constraint equation is 

[j(f + vHR)} 11 ' 1 = -AnGaR. (8) 

Solving Eqs.(4) and (5) we can obtain the time evolution of a and vn in the shell, and v l is derived 
using Eq.(7). These values of a, v l and v 11 satisfy the condition (8). The initial condition is given 
as a form of {H l )i = {H ll )i at an initial epoch (P)i such as the recombination epoch. 

The background models in V 1 and V 11 are rewritten using 

y 3 ee a 3 /(a y, r 3 ee H 3 t 3 , A^ ee \Ac 2 /(H 3 f, nj ee ^-(poY , (9) 

where denotes the present epoch. Eq. (3) is given as 

dy 3 /dr 3 = (y 3 )- 1 / 2 P J (y 3 ), (10) 

where 

PM) ee [nj + \i(y 3 f + (i - nj - W'] 1 ' 2 (11) 

and (ao)- 7 is given by 



(a tf ) J = 1/Vl-^o-^o- (12) 
The conformal times rf are defined by 

rf ee ^Jl-ni-Xij" dy/^P^y)]. (13) 
The solution of Eq. (10) in the case of k 1 = k u = — 1 and A = are expressed as 

n j 

1 - (coshV-1), (14) 
2(1 — SZ J 

- J - ^ (sinh^ - T] 3 ). (15) 



2(1 - nj)3/2 

In the case of nonzero A, Eq. (10) is solved numerically. 

Eqs. (4) and (5) were solved by Sakai et al. (1993) and it was shown in the case of ~ 1 that 
the present value of v 11 at a/cii > 100 is < 100 km/s, as long as the shell starts with the vanishing 
initial velocity (v J )i = 0. 
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Here the initial condition (H L )i = (H lI )i at initial epoch y- 7 = (y^)i is expressed as 



H, 



where we have 



„ (ylr 3/2 [nl + xl {y]f + (1 - n l - XDy]] 1 ' 2 

= H$(tf)-W[t$ + A^) 3 + (1 - O? - A?)yF] 1/2 , 



\H U J 



from Eq. (9). If we eliminate SIq 1 , Aq 1 from Eq. (16) using Eq. (17), we obtain 

[4(1-,,")- $)\i-yiM- + \( Hlh2 

Since y\ ~ y? 1 ~ 10~ 3 (<< 1), we assume 
where e « 1. Then we get from Eq. (18) 



ii 



= 0. 



^ = Pr)'W), 

Po v 



l-(^V 



H \2 



n \3 



(16) 
(17) 

(18) 



(19) 



(20) 



If we give Slg,Ag,e and yj l /yj, therefore, we can obtain H^/Hl for Eq. (20) and derive p$ / p\ 
from Eq. (20) and Vt^/VL\ and X l ^/X\ from Eq. (17). From Eqs. (17), (19) and (20), we find that 

#o7#o < 1 and fi?/ n o > !> if 



1 > > 1 + 

2/1 



2/1 



1 -(!+*)% 

2/1 



4- 



(21) 



In the case A = 0, we have an example for f2 = 0.2 

fi? = 0.56, H$/Hl = 0.80, PqVpo = 1-8, yf/y] = 1.2, e = 4.1. 
In the case A / 0, we have an example 

X\ = 0.672, A r = 0.43, fijj = 0.3, f^ 1 = 0.563, fl^/fTj = 0.80, e = 0.64, 

so that 



n l + A = 0.872, fi^ 1 + Aq 1 = 0.993. 



ii , \H 



(22) 

(23) 
(24) 
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2.2. Bulk motion 

Now let us consider the velocity field around an observer O in V 1 at the point with Iq = 
( a X)o << ( a x)i- Since Iq is much smaller than the curvature radius, we can approximately neglect 
the spatial curvature around him. Then he has the relative velocity 

Av = (Hi - Hl% (25) 

to matter in the outer region in the direction of the X axis. If Hi = lOO/i km/sec/Mpc, Hf = 
0.82HI and / = 40/i" 1 Mpc, we have Av = 720 km/sec. 

Here consider a galaxy G with the radius coordinate x an d angle (p. Then the relative velocity 
of G to matter in the outer region is 

Av G = (Hl-Hl I )(a oX ) (26) 
in the radial direction from the center C of the inner region (cf . Fig. 3) . 

This velocity can be divided into the X component (Avc)x and the line-of-sight component 
(Avg)ls with respect to the observer O as follows, noticing that the angle ZGOX (= 4>) satisfies 
the relation 

sin <p = x sin ip / {x 2 + Xo-2xXo cos ip} 1/2 , cot <p = (cos tp - Xo/x)/ sin ip. (27) 
Their values are 

(Avg)x = Avq sin(0 — (f)/ sin <p = Avq, (28) 

(Av g )ls = Awe sin 92/ sin = Av G [l - 2(xo/x) cos 99 + (xo/x) 2 } 1/2 - &v G . (29) 

That is, the X component (Av G )x is constant and equal to Avo, and another component, which is 
in the line-of-sight direction from the observer, is nearly equal to Av G . Because the first component 
(Av G )x is independent of the position, it can be interpreted as the peculiar velocity v p of galaxies 
which represents the bulk motion: 

v p = (Av G )x = Av . (30) 



3. Redshift formula and the CMB anisotropy 



The wavevector in the inner and outer regions in the plane of = -it/ 2 is obtained by solving 
null-geodesic equation and expressed as 



(k°y = dp/d\ = ^w{, 



(k x ) J = dx 3 /d\ = ± 



3 1 



a J sinh x J 



1/2 



(k*) j = dp/dX = d j /(a j sinh X j ) 2 , 
where j = I and II, and A is an affine parameter. 



(31) 
(32) 
(33) 
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3.1. Light paths 

In the inner region V 1 , it is assumed that at the present epoch (t l = tf,) all rays reach an 
observer at the point O with x = Xo, = ir/2 and ip = in the X axis, and the angle between the 
rays and the X axis is <p. Then we have 



where 



4> = 4>i and 7r — 0i , 



4n = sin- 1 h^irl « vr/2) 
l % smh Xo J 

For <p = <p\ the rays are expressed 



as 



G( X ') S cosh"' (SS*2d) - cosh"' = - 



where 



r/ 1 is defined by Eqs. (13) and (14), and r]\ is equal to if 

hV\t* A\ — tt — /Tt-, TTtro htnro 



1/2 



l v ya )-> / j 
r7g is equal to r/ 1 at present epoch (y 1 = 1). 

For <p = 7r — <pi we have 



where 



G^X 1 ) = — »7o + ^) for Vo>V I > Vm, 

= -if -V0 + 2 Vrn, for if < T) m , 



(34) 
(35) 

(36) 
(37) 



I , -i /coshxo\ 

In the latter case, x 1 has the minimum value (i.e., k x = 0) at if = r] m . 
At the boundary rf = rf and X 1 = Xii therefore, we obtain 

I I , u-i f cosh xl \ , -i /coshx[\ 
% = % ± cosh - cosh ^-^j-j 

for (j) = ( | , respectively. In the outer region V 11 , we have 

\*-<PiJ 

G(x") s cosh- - cosh- (£=■£<?) = „« - 

where r/ n is given by Eqs. (13) and (14), and rff and x! 1 are the values at the shell, 
recombination epoch we have 

i / cosh vll: _ \ i /cnsh v?\ 



(38) 
(39) 



(40) 



rSL = Vf ~ -h- (£2^) + cosh- (^) . 



(41) 
At the 

(42) 
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The junction of wavevectors at the boundary is expressed as 



Eqs. (43) and (45) give 



(a/ 9 (pip k^) = (a/ 3w • 



(43) 

(44) 
(45) 



a o\ r i 



fl 



. (46) 
ai/ Vai/ 

d 1 = d n , (47) 

respectively, where we used the relation R = (a/) 1 = (af) u . The conditions (43) and (45) arc 
evidently consistent with Eq. (44), since k^ is a null vector. 



3.2. Redshift formula 

Now let us derive the redshift formula for rays which are emitted at the recombination epoch. 
Here this epoch is defined as the time of the radiation temperature T r = 10 3 (T r )o in the region V 11 , 
where (T r )o is the present temperature (~ 2.7 K). The total redshift factor (1 + z TCC ) is calculated 
as the product of two redshift factors which are caused in the two regions V 1 and V 11 . 

First we assume that the shell is comoving and later the correction due to the motion of the 
shell is examined. If we consider a virtual observer at the center C (% = 0), a light ray which is 
received by him at if = fj^ is expressed as 

f\l - V 1 = x\ vl - fj\ = X \ (48) 

in V 1 , and 

Vrcc V Xrec X ? Vrec Vl Xrec Xl 

(49) 

in V 11 , when the ray is emitted at the recombination epoch. 
The redshift factors are 



a l {f)\) y 1 (fj[) 

in V 1 , and 



1 + ^ = ^IT = ^K (50) 



1 + zH a u (fj?) y ll {nY) 



1 + z? a^fjllc) y u (fjH c ) 
in V 11 . The junction condition in Eq. (6) gives z\ = z} 1 for v l = v 11 = 0, so that 



(51) 



1 + ^rec — n / -r/? 1 tI -T \ • (52) 



y l \mlc)y l {n[) 
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Here we specify the value of z\ as z\ = 0.067 ~ 0.1. Then r]\ and x\ are determined from Eqs. (48) 
and (50), and fj^ and x? are determined using the relations 

a\y\rj\) sinh X \ = a^HvY) sinh X ? (53) 
and | 

JO JO 

which are obtained from Eqs. (2) and (6). If zll c is moreover specified (in the following we take 
the value zll c = 10 3 — 1), y(fj iec ) or ?7 rcc is determined using Eq. (52). 

Next we consider the observer at O (with \ l = Xo an d <P = 0)- The above determined 
x{ (j = I and II) and r]ll c (= fjH c ) are used also for rays reaching this observer, to identify the 
position of the shell and the recombination epoch. It should be noted that i]H c depends only on 
the temperature at the recombination epoch and is independent of the existence of the inner region 
V 1 . 



In V 7 we have the redshift 



and in V 11 we have 



1 + ^ = ^7TT = ^fTT ( 55 ) 
a l {r)[) y L (r]\) 

1 + zH a 11 ^ 1 ) y 11 ^ 1 ) 



1 + z? a^H) y n (r? n c) - ( 56 ) 

For a given \\-> we obtain r]\ by solving Eq. (40) and obtain z\ from Eq. (55). The junction 
condition in Eq. (6) gives z\ = z} 1 , and is related to r)\ using the relations (53) and (54). Then 
the product of Eqs. (55) and (56) reduces to 

which is given as a function of angle 

Now let us consider the case when the shell is not comoving, i.e., v 1 / 0, v 11 / 0. As was shown 
by Sakai et al., the velocities are < 200 km/sec and (v^ /c) 2 < 10~ 7 for j = I and II. Accordingly, 
we can assume 7 J = 1 for j = I and II, so that the condition z\ = z} 1 and Eq. (54) hold. In order 
to take the shell motion into account, we assume the relation 

X I i = (x I i)^=o + « I [^-(^)0=o] (58) 

in Eq. (40) for arbitrary (j), where (xi)<^=o an d ( 7 ?!)<^=o are their values for the ray incident in the 
direction of <p = 0. Subsequent calculations are same as the calculations for v l = 0. 

The derivations of zll c in the case of double shells and in the inhomogeneous model with an 
intermediate self-similar region are shown in Appendixes A and B. 
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3.3. CMB anisotropy 



The values of zll c are numerically calculated for < (f> < tt, and their dipole and quadrupole 
moments are derived. When z iec (= zll c ((f))) is given, the temperature T(4>) of the cosmic back- 
ground radiation is proportional to 1/(1 + z rcc ), and the dipole moment D and quadrupole moment 
Q are defined as 

f7T / i 2-7T 



D = 



Q = 



f'TT f'ZTT 

/ / (1 + Z^Y, 
JO JO 



io sin 



/*7r f'2n; 

/ / (1 + z xcc )~ x Yio sin 
Jo Jo 



I < (1 + Z r cc) 1 >, 

/ < (1 + z rcc )- L 



>, 



where <> means the average value taken over the whole sky , and 



Y 10 (<f>) = 




coscj) Y 20 ((f)) = 



5 /3 



4tt ~ " uvr/ V 47r \2 

The Doppler velocity Vd corresponding to D is given by 

BiE^jr) 1 / 2 !?]. 



cos 
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(59) 
(60) 

(61) 
(62) 



Assuming Iq (= aoXo) = 40(/i I )~ 1 Mpc mainly, D, Q and Vd were derived for various model param- 
eters. Their values in models with a single shell and double shells are shown in Tables 1 and 2, 
respectively. The values in models with an intermediate self-similar region are shown in Table 3. 
In Table 1 the values in the case with a moving shell are shown. 

It is found from Table 1 that (1) Vd/v p is smaller than 0.17 in all cases, (2) Vd is approximately 
proportional to Iq (as well as v p ), (3) Vd is smaller for smaller (Q l — Qq) and (h l — h 11 ), (4) the 
positive cosmological constant plays a role of increasing Vd, and (5) the influence of the shell motion 
on D, Q and Vd is negligibly small. 

If we compare two lines in Table 2 with the corresponding ones (the first line and 8th line) 
in Table 1, the results in the models with a single shell and double shells are found to be quite 
consistent. Moreover, if we compare four lines in Table 3 with the corresponding ones (the second 
line, 4th and 6th line) in Table 1, the results in the models with a single shell and a self-similar 
region are found to be similarly consistent. Accordingly Vd/v p is w 0.1 in all models we treated 
here. 



4. Naive derivation of redshift factors 

For the dipole anisotropy the maximum difference of temperatures and redshifts can be seen 
in a direction (0 = <j>{) and the inverse direction (</> = tt — (f>\). Here let us compare the redshifts 
in the directions = and tt appearing in the model with a single shell. The spatial curvature 
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is neglected for simplicity. In the X axis we consider six points 0, A, B, C, D, and E, for which 
X = Xo, Xi, Xi + 2^0, 0, -xi, and ~xi + 2xo, as shown in Fig. 4. 

Points B and D have the equal distance from the observer's point O. 

The redshifts z rc c(0) and z TCC (ir) of the rays from another points P and P' at the recombination 
epoch to the observer at O in the directions <j) = and ir are divided into three steps 

(P -> B, B -> A, A -> O) and (P' -> D, D -> E, E -> O), 

respectively. That is, 

1 + ^rcc(O) = (1 + 2fp B )(l + ZBa)(1 + Z AO ), 1 + Z tcc (tt) = (1 + *P' D )(1 + *De)(1 + *Eo)- (63) 

Among these three steps, the first and third steps have equal redshifts evidently : zpb = zp'D, zao = 
zeo- in the processes (B — > A) and (D — ► E), we have the redshifts due to the cosmic expansion, 
the Doppler shift and the gravitational redshift. The expansion redshift factors in V 11 and V 1 are 

[l + <(2a oX o/c)], [l + fl3(2aoxo/c)], (64) 

respectively. The Doppler shifts at A and D due to the relative velocity between V 1 and V 11 frames 
are 

[1 --(Hi- <)a (xi - Xo)], [1 - -(H l Q - <)a (xi + Xo)], (65) 
respectively. The gravitational redshifts are represented by potentials f/'BA and V'de given by 



■0BA = 



_ 4^G(a ) 3 rlPo(xi) 3 , P I (xi+Xo) 2 Xo 



c 2 



+ 



3 Xi + Xo Xi + Xo 



_ 47rG(a ) 3 rl Po(xi) 3 Po(Xi ~ Xo) 2 Xo 
Vbe = 5— 



3 Xi + Xo Xi - Xo 



(66) 
(67) 



Accordingly the total redshift factors are 



1 + zba = [l + <(2a Xo/c)][l-^(^-0«o(xi-Xo)][l+^BA] 

= l-f(^-Oxi + |(^ + <)X0+^BA 

-(^y H l\Hl - <)xo(xi - Xo) + • • • , (68) 

l + z DE = [l + Hl(2aoXo/c)][l-±(Hl-H$)ao(xi+Xo)][l+ifoE] 
= i_^ (jff i_ F n) Xl+ ^ ( ^i + jff n )xo + V)DE 

- (^^M " <)xo(xi + Xo) + • • • • (69) 
The difference of these factors reduces to 

Z mc (0) - Z tcc (tt) ~ Z B A - ^DE 
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V'ba - V>de + (^) "(Hi - Oxo [(Hi + Oxo + (ffS - OxiJ • (70) 



That is, the main terms in (1 + zba) and (1 + zde) cancel. The ratio 1Z of [z rec (0) — z rcc (ir)\ to the 
relative velocity between V 1 and V 11 is 



n 



(0) - Z rcc (ir) 

(Hi - Hl l )a 0X0 /c 



4W) 2 (Xi + Xo) - tf (Hl) 2 ( X i ~ Xo)] - Hi 1 ) 
-2(Hl + <)xo/c + 2(i*S - Oxi/c, 



(71) 



where we used Eq.(9) for Pq. This ratio is the counterpart of the ratio of the velocity Vd (corre- 
sponding to the dipole moment) to the relative (peculiar) velocity v p , and these two ratios have 
comparable values. 



5. Consistency with the observed large-scale bulk flows 

As was shown in B of §2, the bulk velocities in all positions within the region I are equal in the 
present models, so that the relative velocity (t>Lc) of LG to the cluster frame is only the peculiar 
velocity of LG caused by GA and the nearby superclusters. Moreover, as was shown in C of §3 the 
dipole velocity Vd corresponding to the bulk velocity v p is ~ 0.1v p , and so the total dipole velocity 
v td of LG to CMB is v td = v LG + v d (~ 0.1v p ) = v LG . 

These conclusions are consistent as follows with (1) the observed velocities (vlg) of LG with 
respect to the cluster frames by Giovanelli et al. (1998), Dale et al. (1999) and Riess et al. (1997), 
which are nearly equal to the total dipole velocity (v t d) with respect to CMB, and (2) the bulk 
velocities (v p ) in the SMAC observation (Hudson et al. (1999)) and the LP 10k observation (Willick 
(1999)) are ~ 700 km sec -1 in the nearly same directions: 

According to Dale et al. (1999), we have v LG = 565±113 km sec" 1 and (I, b) = (267°, 26°)±10°. 
On the other hand, v td = 627±22 km sec" 1 and (I, b) = (276°, 30°) (Kogut et al. (1993)). Since both 
directions are nearly equal, the velocity difference is about 60 km sec -1 . This value is comparable 
with Vd (~ 0.1-Up), in the consistent manner with the result in the present models. 

However, they are inconsistent with the observations of Lauer and Postman (1994) and Colless 
(1995) in which vtd and vlg are in quite different directions. At present these observations seem to 
have been ruled out (cf. Proceedings of the International Conference on Cosmic Flows 1999). 



6. Concluding remarks 

In this paper we considered the behaviors of galaxies and light rays in spherically symmetric 
inhomogeneous models consisting of inner and outer homogeneous regions V 1 and V 11 with S7q and 
SIq 1 (> S7g) and Hi and H^ (< Hi), respectively, connected by a single shell. In Appendixes we 
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treated also the models with double shells and an intermediate self-similar region. It was shown 
as the result that when we observe the motion of galaxies at the point O deviated from the center 
C in V 1 , a constant peculiar velocity component v p appears in the direction from the center to 
the observer (C — ► O). Moreover it was shown that the velocity Vd corresponding to the dipole 
anisotropy of CMB radiation is by a factor « 10 small compared with v p . This result may fit for 
the observed situation of the cosmic flow of cluster galaxies, when the scale of the inner region 
and the distance CO are about 200(/i 1 ) -1 and 40(/i 1 ) -1 Mpc, respectively, and when (f^ 1 — Qq) is 
~ 0.5 and (Hq — H^) is ~ 18h l km/sec/Mpc (or h ll /h l ~ 0.82), respectively. This difference 
of the Hubble parameters may be consistent with their recent values due to nearby and remote 
observations, because short and long distance scales give the Hubble constant of Hq > 70 and 
Ho ~ 55, respectively (Branch 1998; Freedman 1997; Sandage and Tammann 1997; Blandford and 
Kundic 1997). A model with multi-shells in the region of 100 ~ 300 Mpc in which the parameters 
change stepwise will be better to reproduce the observed distribution of the Hubble constant. 

It was shown that the present models are consistent with the current observations of large- 
scale flows (Giovanelli et al. (1998), Dale et al. (1999), Riess et al. (1997), Hudson et al. (1999) and 
Willick (1999)), but inconsistent with other observations (Lauer and Postman (1994) and Colless 
(1995)) which may have been ruled out. 

It is interesting and important to study the influences of the above inhomogeneity on the 
cosmological observations such as the magnitude-redshift relation of SNla, the number count of 
galaxies, the time delay for lensed QSOs, and so on. They will be quantitatively analyzed and 
shown in near future. 

In this paper the motion of LG due to GA and superclusters in similar distances was not 
treated, while spherical matter distributions on such scales were analyzed by Humphreys et al. 
(1997). But the approximation of spherical symmetry may not be good to their small-scale matter 
distribution. 



Models with double shells 



The spacetime is divided into three homogeneous regions V 1 , V 11 , and V m , as shown in Fig. 
5. The cosmological parameters in V- 7 are shown as Q J , \ J and so on, where j = I, II, III. 

Eqs. (1) - (3) and (9) - (15) hold also in region III as well as in regions I and II. The junction 
condition at the second boundary between V 11 and V m have similar expressions to those at the 
first boundary: 

(Al) 

iiiji 

(A2) 



7 III d(47r J RV)/dt ili = AttR^vp 



R^a ll f u 



in 



a m f Ul , 



}2 2„ 



d( 7 m v m )/dt m = - 7 m v m H m + 2vrGa - h 2 v 2 p] m /a, 



(A3) 
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dt u /dt U1 = 7 n /7 in , 
7 ii (/ /ii v ii + H n R) = 7 iii (/ /iii v iii + H m R)j 



and 



[ 7 (/' + v#-R)] ni ' n = -AirGaR, 

where [$] in - n = $ m - $ n . 

Similarly to Eqs. (19) and (20), moreover, we have 

-Hi /? ,m x3 
Po^ 1 



(A4) 
(A5) 

(A6) 



V Hi 1 ) \ v 11 J 



Po v iff y 



U n / J w 11 / 



+ A, 



n 



1 



'0 > 



(A7) 
(A8) 



where e 



As for light rays, the equations in V 1 and V 11 are same as those in Sec. 3, and equations in 
V m are common with those in V 11 . At the second boundary (r? = 772, X = X2) we have 



/a \ n n /«o\ m in 



d 



ai- 

• u -d n \ 



(A9) 
(A10) 



where a2 = 0(772)- 

Under the assumption that the shells are comoving (i.e., v 1 = v 11 = v 111 = 0), we obtain the 
following redshift formulas: 

(AH) 



i + *i = £(rf) 



and in V 11 we have 



1 + 4 1 
1 + ^ 
in 



y u (v¥Y 



1 + 4, 
1 + 4 11 



(A12) 

y^m- (A13) 

From Eqs. (6) and (A4) we have two equations z\ = z} 1 and = z^ 1 . Moreover, (^jXi 1 ) arc 
related to (rj\, \\) using Eqs. (53) and (54), and (f?^ 11 , x^ 1 ) are related to (r/J, 1 , x^) using the similar 
equations 

a$y l \r$) sinh^ 1 = a l ^y m (^ 1 ) sinh X2 n (A14) 



and 



JO JO 



(A15) 
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First we consider the virtual observer in the center (x = 0). Then we have 

1 + ^ = — ^- (A16) 

l + ^II yII fe )H' ^ 



i + 4 n y m (^ec) m ' (A18) 

where ffo — fj\ = \\ an d V2 ~ = x\ l ~ X?- Accordingly, we can determine x^x!; 1 an< ^ ^rcc by 
specifying his redshifts z\, z^ and (= 10 3 — 1), as in Sec. 3. 

For an observer at O (with r] = t]q and x = Xo)j we obtain 



y ii {r] u )y m iv m 



x ^ ^rcc t/ T\...TT/_TT\..TTT/_TTT\ ' 



where r/{ , 772 , r/2 an d 7/rcc depend on the angle 4>, contrary to those in Eq. (A18). When we fix 
z\, and z^} c and determine x\iX\ l an d Vreci the final redshift z^} c can be obtained as a function 
of angle cp for 77™ = 77™. ^ is to be noted that 77™ is independent of the existence of the two 
shells. 



B. Models with an intermediate self-similar region 

The line-element is expressed in the form (Tomita 1995, 1996) 

ds 2 = _ c 2 dt 2 + S 2 {t r J (1 + rg V g ) 2 dr 2 +r 2 dn 2) (m) 

I 1 — ka(r)r z j 

where 

a(r)/oo = V(n) 2 , 1/r 2 , and l/(r 2 ) 2 (B2) 

for the inner homogeneous region V 1 , the self-similar region V 11 , and the outer homogeneous region 
V m , respectively, which are shown in Fig. 7. Assuming the matter pressure- free and comoving, the 
solutions are described using the Tolman solution. The two boundaries r = n and 7*2 are exactly 
comoving. The scale factor S in the case A = is given by 

S/So = ^_ (cosh 7,-1), (B3) 

and 

n = (cH 1 /S )\ H =[(dS/dt)/S} . (B4) 
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The time variable t is in the inner, self-similar, and outer regions expressed as 

^0 



2(1 - ^o) 3 /2 



(sinhr/ — rj), 



ct cHo-^n)- 1 ^ 



r 2(1 - O ) 3 / 2 



(sinhry — rj), 



H t 



T2 



n [2(1 - f) )3/2] 

respectively.he constant «o and n are given by 



(sinhr/ — rj). 



yi - n 
n 2 (i + zi) 



i-n + ^(1 + ^) + 



1 + n 



QZl 



(B5) 

(B6) 
(B7) 

(B8) 



and 



r 1 = [(n ) 1 ' 3 /(l-n ) 1/2 WaS- (B9) 

If we define the local Hubble parameter if (in the transverse direction) in the t = to hypersurface 

as 

H = (S/S) t=t0 = [(l/mS^/S] t=t0 , (BIO) 
H in V 1 is equal to Ho and we have in V 11 and V 1 

sinhr/(sinhry — fj) 



Hto 



(cosh rj — l) 2 



where rj = rj t=to satisfies 



#0*0 — ( — , — ) 



Or 



n'nJ [2(1 -O ) 3 /2] 



(sinhr/ — 7/). 



(Bll) 



(B12) 



for (V , V , respectively. Present matter densities in the three regions are 



3Ho 2 Qq 



8vrG 



1, 



1 (t\ 



3 V r 



1 - fS^/S 



-1 / \ 2 
ri 

^2 



(B13) 



for j = [I, II, III], respectively, and SIqj is defined as 



Qoj = Poj(to)/ 



3H 2 
8vrG 



(B14) 



The behavior of light rays and redshift formulas were derived in the previous paper(Tomita 
1996). If we specify the redshifts z±, zi and z rec (= 10 3 — 1), the radii r\ and r2 of two boundaries 
and the value of r] TCC (= fj TCC ) are determined, and we can calculate z rec as a function of (f>, which is 
the angle between the incident direction and the X axis. 
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In the following we present additional explanations and corrections for the formulas in the 
Appendix A of the previous paper(Tomita 1996) : In V 1 and V 111 we have 



¥ = ±(S /S 2 ){[1 + ao(r/ ri ) 2 ][( Wl f - d 2 /(-V) 2 ]} 1 ^ , 

V = -(S /S 2 ){[1 + a (r/r 2 ) 2 ][(w 2 ) 2 - d 2 /(-V) 2 ]} 1/2 , 
respectively. The second ray in V 1 has 



hi = cosh 1 ■ 



l + (vWVn) 



1/2 



/ho > + {(v ~ Vo), ~(V ~ 2r lm + Vo)} 



for {rj > r] m ,rj < rj m }, respectively. In this paper the incoming angle is 



sin 



Auxiliary functions V 111 are expressed as 



M(0 



N(0 



(M ) 2 -2d 2 J^ dZ 



MO 



1/2 



M 2(1 -fi ) cosh??-l /M 



d 



and 



c 



fio^o 2(1 — cosh?]) + ?7sinh?7 V d 
oq 2(1 — cosh rj) + r/sinhr/ 



2V 



i/T+~ao (coshr? — l)(sinh?7 — rj) 



(B15) 
(B16) 

(B17) 



(B18) 



(B19) 

(B20) 
(B21) 

(B22) 



The constants w\ and w 2 are connected using the junction conditions (fc ) 1 = (A; ) 11 at r = r± 
and (k°) u = (k°) m at r = r 2 . 

The author thanks N. Sugiyama and N. Sakai for helpful discussions, and a referee for kind 
comments about past related works. This work was supported by Grant-in Aid for Scientific 
Research (No. 10640266) from the Ministry of Education, Science, Sports and Culture, Japan. 



REFERENCES 

Blandford, R. D. and Kundic, T. 1997, The Extragalactic Distance Scale, M. Livio, M. Donahue 
and N. Panagia, London: Cambridge University Press, 60 

Branch, D. 1998, ARA&A, 36, 17 

Broadhurst, T. J., Ellis, R. S., Koo, D.C., Szalay, A. S. 1990, Nature, 343, 726 



-18- 



Colless, M. 1995, AJ, 109, 1937 

Dale, D. A., Giovanelli, R., and Haynes, M.P. 1999, ApJ, 510, Lll 

Dressier, A., Faber, S.M., Burstein, D., Davies, R. L., Lynden-Bell, D., Terlevich, R. J., and 
Wegner, G. 1988, ApJ, 313, L37 

Einasto, J. et al. 1997, Nature, 385, 139 

Freedman, W. L. 1997, Critical Dialogues in Cosmology, N. Turok, Singapore: World Scientific, 92 

Giovanelli, R., Haynes, M. P., Freudling, W., da Costa, L. N., Salzer, J. J., and Wegner, G. 1998, 
ApJ, 505, L91 

Hudson, M. J., Smith, R. J., Lucey, J. R., Schlegel D. J. and Davies, R. L. ApJ, L79; Hudson, M. 
J., Lucey, J. R., Smith, R. J., and Steel, J. 1997, MNRAS, 291, 488 

Humphreys, N. P. et al. 1997, ApJ, 477, 47 

Israel, W. 1966, Nuovo Cim. 44b, 1 

Kogut, A. et al. 1993, ApJ, 419, 1 

Landy, S. D., Shectman, S. A., Lin, H., Kirshner, R. P., Oemler, A., and Tucker, D. 1996, ApJ, 
456, LI 

Lauer, T. R. and Postman, M. 1994, ApJ, 425, 418; Postman, M. and Lauer, T. R. 1995, ApJ, 440, 
28 

Lynden-Bell, D., Faber, S. M., Burstein, D., Davies, R. L., Dressier, A., Terlevich, R. J. and 
Wegner, G. 1988, ApJ, 326, 19 

Maeda, K. 1986, Gen. Rel. Grav., 18, 9 

Nakao, K., Gouda, N., Chiba, T., Ikeuchi, S., Nakamura, T. and Shibata, M. 1995, ApJ, 453, 541 
Riess, A. G., Davis, M., Baker, J., and Kirshner, R. P. 1997, ApJ, 488, LI 
Sakai, N., Maeda, K. and Sato, H. 1993, Prog. Theor. Phys., 89, 1193 

Sandage, A. and Tammann, G.A. 1997, Critical Dialogues in Cosmology, N. Turok, Singapore: 
World Scientific, 130 

Tomita, K. 1996, ApJ, 461, 507. 

Tomita, K. 1995, ApJ, 451, 541; Tomita, K. 1996, ApJ, 464, 1054 for Erratum 
Willick, J. A. 1999, ApJ, 522, 647 



-19- 

Zehavi, L, Riess, A. G., Kirshner, P. and Dekel, A. 1998, ApJ, 503, 483 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-20- 



Fig. 1. — A model with a single shell, z and z are the redshifts for observers at O and C. 

Fig. 2. — A schematic diagram of Hubble and density parameters. A solid and dotted lines denote 
Hq and f^o, respectively. 

Fig. 3. — Components of the bulk velocity. 

Fig. 4. — A diagram for the naive derivation. 

Fig. 5. — A model with double shells. 

Fig. 6. — A schematic diagram of Hubble and density parameters in a model with double shells. A 
solid and dotted lines denote Hq and CIq, respectively. 

Fig. 7. — A model with the self-similar region. 
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Fig. 8. — A schematic diagram of Hubble and density parameters in a model with the self-similar 
region. A solid and dotted lines denote Hq and Qq, respectively. 
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Table 1: Dipole and quadrupole moments and the velocity Vd in the single-shell models. 









X 11 


h 1 


h ll /h l 


D (xlO 4 ) 


Q (xlO 5 ) 


Vd (km/si 


0.2 


0.56 


0.0 


0.0 


0.7 


0.82 


5.56 


-2.87 


81.8 1 


0.2 


0.88 


0.0 


0.0 


0.7 


0.82 


8.27 


-4.27 


121.5 1 


0.2 


0.56 


0.0 


0.0 


0.7 


0.90 


5.38 


-2.79 


79.1 1 


0.2 


0.88 


0.0 


0.0 


0.7 


0.90 


8.21 


-4.23 


120. 7 1 


0.3 


0.56 


0.0 


0.0 


0.7 


0.82 


3.90 


-2.04 


57.3 1 


0.3 


0.88 


0.0 


0.0 


0.7 


0.82 


6.66 


-3.39 


97.9 1 


0.2 


0.56 


0.672 


0.430 


0.7 


0.82 


7.92 


-3.74 


111.4 1 


0.2 


0.56 


0.0 


0.0 


0.7 


0.82 


5.55 


-1.91 


81.6 2 


0.2 


0.88 


0.0 


0.0 


0.7 


0.82 


8.21 


-2.68 


120. 7 2 


0.2 


0.56 


0.0 


0.0 


0.7 


0.82 


8.07 


-6.35 


118.6 3 


0.2 


0.88 


0.0 


0.0 


0.7 


0.82 


12.31 


-9.62 


181.0 3 


0.2 


0.56 


0.0 


0.0 


0.7 


0.82 


5.57 


-2.87 


81.8 4 


0.2 


0.88 


0.0 


0.0 


0.7 


0.82 


8.27 


-4.27 


121.6 4 



1 z\ = 0.067, l = m^h 1 )- 1 . 
2 z\ = 0.1, l = 40(h 1 )- 1 . 
z z\ = 0.067, l = 60(/i 1 ) -1 . 

A z\ = 0.067, l = 40(h l )-\ and the shell velocity is 200 km/sec. 
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Table 2: Dipole and quadrupole moments and the velocity Vd in the double-shell models in the case 
of X J = and l = 40(/i 1 )- 1 . 

ttl tig n^ 1 h l h ll /h l h m /h l D (xlO 4 ) Q (xlO 5 ) v d (km/sec) 
~~02 036 056 7 092 082 5.55 -2.37 81.6 1 

0.2 0.36 0.56 0.7 0.92 0.82 5.52 -1.44 81. 1 2 



1 z\ = 0.067, zf = 0.1. 
2 z\ = 0.1, zf = 0.2. 
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Table 3: Dipole and quadrupole moments and the velocity Vd in the model with a self-similar region 



n the 


case 


of\i 


= and Iq 


= 40(/i 1 )- 1 . 






nl 




h l 


h m /h l 


D (xlO 4 ) 


Q (xlO 5 ) 


Vd (km/sec) 


0.2 


0.89 


0.7 


0.81 


8.40 


-4.41 


123.4 1 


0.2 


0.72 


0.7 


0.83 


7.40 


-2.61 


108.8 2 


0.3 


0.91 


0.7 


0.84 


8.08 


-4.11 


118.8 1 


0.3 


0.80 


0.7 


0.86 


6.94 


-2.25 


102.0 2 



x z\ = 0.067, zl 1 = 0.1. 
2 z\ = 0.1, zf = 0.2. 



